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I. INTRODUCTION 

In recent years many Collaborations have presented 
very interesting results of neutrino oscillations that aim 
at a precise determination of mass-squared splitting and 
mixing angles in disappearance and v e appearance 
measurements. Since modern experiments, such as Mini¬ 
BooNE, SciBooNE, ArgoNeuT, MINERVA, and T2K, 
are performed with detectors made of medium-heavy 
nuclear targets, e.g, Carbon, Oxygen, Argon, or Iron, 
a clear understanding of neutrino-nucleus interactions, 
where all nuclear effects are well under control, is required 
for a proper analysis of experimental data. The recent 
progress, the questions and challenges in the physics of 
neutrino cross sections are reviewed in [IJ-Q. 

The first measurements of the charged-current 
quasielastic (CCQE) flux-averaged double-differential 
cross section on 12 C in the few GeV region by 
the MiniBooNE Collaboration 0,0 have raised exten¬ 
sive discussions that effects beyond the impulse approx¬ 
imation (IA) may play a significant role in this energy 
domain |6h17I|. Multinucleon mechanisms, 2p-2h excita¬ 
tions and meson-exchange currents (MEC), and also RPA 
corrections appear essential to describe the data. Models 
based on the IA, which make use either of a realistic spec¬ 
tral function obtained within a nonrelativistic framework 
Em or of a relativistic IA (RIA) [2(1 E l , generally un¬ 
derestimate the MiniBooNE CCQE cross sections. Only 
the relativistic Green’s function (RGF) model is able to 
give a good description of the data [24j. Although un¬ 
der many aspects based on the RIA, the RGF model can 
recover contributions of final-state channels that are not 
included in other models based on the RIA. These con¬ 
tributions are recovered by the imaginary part of the rel¬ 
ativistic optical potential that is used in the RGF model 
to describe final-state interactions (FSI). 

The RGF model was originally developed within a non¬ 
relativistic [25|, E| and then a relativistic [E, E| frame¬ 
work to describe FSI in the inclusive quasielastic (QE) 
electron scattering. The model was successfully tested 
against electron scattering data [E . E, EHUI and it 
was later extended to neutrino-nucleus scattering, both in 


the charged-current (CC) (E, E - El and in the neutral- 
current (NC) [E E | sector. Although different, the two 
situations present many similar aspects and the exten¬ 
sion to neutrino scattering of the electron scattering for¬ 
malism is straightforward. In the QE kinematic region, 
where the nuclear response to an electroweak probe is 
dominated by single-nucleon scattering with direct one- 
nucleon emission, a reliable description of the FSI effects 
between the ejected nucleon and the residual nucleus is 
very important for the comparison with data. In the 
RGF model FSI are described in the inclusive scattering 
consistently with the exclusive scattering by the same 
complex optical potential (OP) and the components of 
the nuclear response are obtained in terms of matrix el¬ 
ements of the same type as the distorted wave impulse 
approximation ones of the exclusive (e, e'p) process, but 
involve eigenfunctions of the OP and of its Hermitian 
conjugate, where the opposite sign of the imaginary part 
gives in one case an absorption and in the other case a 
gain of strength. In the exclusive scattering, where only 
one channel is considered, the imaginary part gives an ab¬ 
sorption that accounts for the flux lost to other channels. 
In the inclusive scattering, where all elastic and inelastic 
channels are included, the imaginary part redistributes 
the flux in all the channels and in the sum over all the 
channels the total flux is conserved. More details on the 
model can be found in our previous papers (EIE . |32j. 

In other approaches based on the RIA, FSI are in¬ 
cluded in the emitted nucleon state with real poten¬ 
tials, either retaining only the real part of the relativistic 
energy-dependent complex OP, or using distorted waves 
obtained with the same relativistic energy-independent 
mean-field potential considered in describing the initial 
nucleon state (RMF) E E 0. In the relativistic 
plane-wave impulse approximation (RPWIA) FSI are ne¬ 
glected. 

The results of these different descriptions of FSI have 
been compared in |30:| for the inclusive QE electron scat¬ 
tering, in g3| for the CCQE neutrino scattering, and 
in Hi E| with the CCQE and NC elastic MiniBooNE 
data. Both RGF and RMF models can describe success¬ 
fully electron scattering data and their related scaling 
functions. Both models are able to provide a satisfactory 
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description of the CCQE MINERVA data [H, |4lj] . In the 
case of the MiniBooNE CCQE data, both models repro¬ 
duce the shape of the experimental cross sections, but 
only the RGF gives cross sections of the same magnitude 
as the experimental ones without the need to increase 
the world average value of the axial mass Ma [mm. 
The larger RGF cross sections are due to the translation 
to the inclusive strength of the overall effect of inelastic 
channels, including rescattering and some multinucleon 
contributions, that are recovered in the model by the 
imaginary part of the relativistic OP and that are not 
included in the RMF and in other models based on the 
IA. 

The optical potential is a powerful tool to recover and 
include important contributions. The availability of phe¬ 
nomenological relativistic OP’s, obtained through a fit 
of elastic proton-nucleus scattering data, is essential to 
make RGF calculations feasible, but the use of a phe¬ 
nomenological OP does not allow us to disentangle the 
role of a specific inelastic channel and can therefore intro¬ 
duce uncertainties and ambiguities in the interpretation 
of the RGF results. The imaginary part can recover, 
to some extent, contributions beyond direct one-nucleon 
emission, such as, for instance, rescattering of the outgo¬ 
ing nucleon and some multinucleon processes, which can 
be included in CCQE measurements, but the RGF model 
is based on the use of a one-body nuclear current and does 
not contain MEC mechanisms that in other models have 
been found to be significant. On the other hand, the 
imaginary part of the OP can include pion-absorption 
and pion-emission processes, that should have already 
been subtracted in the MiniBooNE analysis. It has been 
written in [3] that the good agreement of the RGF re¬ 
sults with the MiniBooNE data “should be interpreted 
with care” and that “it would be very interesting to con¬ 
front the RGF results with the fully CC-inclusive data” 
Q. The comparison with the fully CC-inclusive data, 
which include also pion production, is the motivation of 
the present paper. 

The T2K collaboration has recently published [42] new 
results on the CC-inclusive double differential cross sec¬ 
tion on 12 C, which includes also pion production, and of 
the CCQE cross section [43[. The T2K energy range 
is the same as for MiniBooNE, the beam peaks at ~ 600 
MeV, similar to that of MiniBooNE, but it is significantly 
narrower and receives almost neglible contributions for 
energies larger than 1 GeV. In view of these differences, 
the analysis of T2K data is another useful and indepen¬ 
dent test for a theoretical description. 

In this paper we compare the results of the RGF model 
with the CC-inclusive v ft and v e T2K cross sections on 
12 C EHH. As a first step, we compare our results with 
the CCQE i/fj, — 12 C cross sections measured by the Mini¬ 
BooNE and T2K collaborations, for which there is con¬ 
sistency between the two experiments within the current 
statistical and systematic uncertainties [43|. Then, we 
consider the flux-averaged CC-inclusive u lx and v e differ¬ 
ential cross sections from T2K with the aim to investigate 
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Figure 1. Total CCQE v ^.- 12 C cross sections per target 
neutron versus the neutrino energy. The experimental data 
are from MiniBooNE Q and T2K [43| . 


whether the effects of the inelastic channels recovered in 
the RGF model by the relativistic OP give or do not give 
enough strength to reproduce these data. 
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II. RESULTS 

In all the calculations presented in this work we have 
adopted the standard value for the nucleon axial mass 
Ma = 1-03 GeV/c 2 . The bound nucleon states are 
taken as self-consistent Dirac-Hartree solutions derived 
within a relativistic mean field approach using a La- 
grangian containing a, w, and p mesons |46|-|50j]. We 
have used two different parametrizations for the rela¬ 
tivistic OP of 12 C that is adopted in our RGF cal¬ 
culations: the Energy-Dependent and A-Independent 
EDAI (where the E represents the energy and the A 
the atomic number) OP of (HTJ , and the more recent 
Democratic (DEM) phenomenological OP of [52;]. The 
EDAI OP is a single-nucleus parametrization, which is 
constructed to better reproduce the elastic proton- 12 C 
phenomenology, whereas the DEM parametrization is a 
global parametrization, which depends on the atomic 
number A and is obtained through a fit to more than 
200 data sets of elastic proton-nucleus scattering data on 
a wide range of nuclei and that is not limited to doubly 
closed shell nuclei. In comparison with electron scatter¬ 
ing data, the DEM parametrization produces in general 
good results for doubly magic nuclei and less good but 
still acceptable results for nuclei with a number of nucle¬ 
ons far from the magic numbers (3lL |53| . 

In Fig. [I] we show our calculated CCQE i/^- 12 C cross 
sections per target neutron as a function of the neutrino 
energy compared with the MiniBooNE Q and T2K data 
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Figure 2. Flux-averaged CC-inclusive double differential z/ M - 12 C cross sections per target nucleon as a function of the muon 
momentum. The data are from T2K [42]. 


{43| . Although the results of these two independent mea¬ 
surements can be consistently compared in the entire 
range of energies, with the only exception of the T2K da¬ 
tum in the energy bin 1 — 1.5 GeV, we observe that the 
average magnitude of the MiniBooNE dataset is larger 
than that of the T2K one. The differences between the 
RGF-EDAI and RGF-DEM results are sizable. These dif¬ 
ferences are due to the different imaginary parts of the 
two OPs, particularly for the energies considered in kine¬ 
matics with the lowest scattering angles and the largest 
kinetic energies of the muon (24J. The RGF-EDAI cross 
section is larger than the RGF-DEM one, in better agree¬ 
ment with the MiniBooNE data and in agreement with 
both MiniBooNE and T2K cross sections within the error 
bars in the entire energy range of the data. The RGF- 
DEM cross section underpredicts the MiniBooNE data 
at low E v and it is in better agreement with the T2K 
data. The RPWIA cross section, which is also shown in 
the figure for a comparison, is similar to the RGF-DEM 
one. We note that other models based on the IA give 


in general results somewhat lower than the RPWIA one 
and therefore lower than the data. 

In Fig. [2] we present the CC-inclusive double differen¬ 
tial z/ m - 12 C cross section d 2 a/(dP^d cos as a function 
of the outgoing muon momentum transfer for four dif¬ 
ferent bins in the scattering angle. The calculated cross 
sections are flux-averaged over the T2K flux [54] and 
compared with the experimental data of (42| . 

The RPWIA results in Fig. [2] are approximately 50% 
lower than the data. Also the RGF results underestimate 
the data. Both RGF-EDAI and RGF-DEM cross sections 
are generally lower than the data, although within the 
error bars for low values of and large angular bins. 
A satisfactory ag reement with the data is obtained with 
the model of [551 ]. which includes np-nh excitations and 
single-pion production. In the RGF model the imaginary 
part of the OP can include the excitation of multinucleon 
channels. We cannot exclude that it can contain some 
contribution due to pion emission, we cannot disentangle 
and evaluate the relevance of this contribution, but in 
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Figure 3. Flux-averaged CC-inclusive u e - 12 C differential 
cross section per nucleon as a function of cos#. Only elec¬ 
trons corresponding to p > 0.550 GeV/c and cos# > 0.72 are 
considered [45| . The T2K data can be found at [44t | . 



Figure 4. Flux-averaged CC-inclusive u e - 12 C differential 
cross section per target nucleon as a function of the elec¬ 
tron momentum. Only electrons corresponding to p > 0.550 
GeV/c and cos# > 0.72 are considered O . The T2K data 
can be found at [44| . 


any case the results in Fig. [2] indicate that this is not 
enough to reproduce CC-inclusive data. 

We note that for the most forward angular bin (0.94 < 
cos < 1.00) the RGF results are significantly smaller 
than the RPWIA ones. In this small bin the trans¬ 
verse and charge-isovector contributions are suppressed 
and the longitudinal response gives the main contribu¬ 
tion to the cross section. In addition, it has been shown 
in f56] that models based on quasi-free scattering cannot 
describe properly this kinematic situation where ~ 1/2 


of the total cross section arises from excitation energies 
below ~ 50 MeV. If we consider that in the RGF calcu¬ 
lations collective effects are neglected, it is not surprising 
that our results are significantly lower than the data in 
the forward angular bin. 

In Fig. [3] and Fig. []] the calculated CC-inclusive v e - 
12 C differential cross sections are displayed as a function 
of the electron momentum and scattering angle, respec¬ 
tively, and compared with the T2K data of |45| . For sake 
of simplicity the calculations have been performed only 
for the reduced phase-space (momentum > 550 MeV/c 
and cost? > 0.72) and not for the full phase-space of 
T2K. The T2K v e beam peaks at ~ 500 MeV and, in 
contrast to the one, extends to energies larger than 1 
GeV. Thus the CC v e cross section at T2K may receive 
also higher-order contributions like two-pion production. 
Also in this case our calculations are significantly lower 
than the data. The RPWIA and the two RGF results are 
very similar in Fig. 0 where the three curves practically 
overlap, while in Fig. d]the differences are small but vis¬ 
ible. The fact that collective effects are not included in 
our model is one of the reasons of the large underesti¬ 
mation of the experimental data at smaller angle in Fig. 
0 


III. CONCLUSIONS 

In this paper we have compared the predictions of the 
RGF model with the CCQE and CC-inclusive and v e 
scattering T2K data. The RGF model is able to give 
a satisfactory description of inclusive QE electron scat¬ 
tering cross sections and of the CCQE MiniBooNE and 
MINERVA data, both for and scattering, without 
the need to increase the standard value of the axial mass. 

The RGF results are usually larger than the results of 
other models based on the IA. In the RGF model FSI are 
described using a complex energy-dependent relativistic 
OP whose imaginary part includes the overall effect of 
the inelastic channels, which give different contributions 
at different energies, and makes the RGF results sensi¬ 
tive to the kinematic conditions of the calculations. With 
a complex OP the model can include all the available 
final-state channels and not only direct one-nucleon emis¬ 
sion processes. The important role of contributions other 
than direct one-nucleon emission has been confirmed by 
different and independent models in the case of CCQE 
MiniBooNE cross sections, but the same conclusion is 
doubtful in the case of MINERVA data. 

The RGF model can include contributions of final-state 
channels like, e.g., rescattering processes of the nucleon 
in its way out of the nucleus, non-nucleonic A excita¬ 
tions, which may arise during nucleon propagation, with 
or without real-pion production, and also some multin¬ 
ucleon processes. These contributions are not incorpo¬ 
rated explicitely in the model with a microscopic cal¬ 
culation: they can be recovered, to some extent, at a 
phenomenological level by the imaginary part of the phe- 
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nomenological OP which is adopted in the RGF calcula¬ 
tion. The use of a phenomenological OP does not allow us 
to disentangle and evaluate the role of a specific reaction 
process. Different available parametrizations of the phe¬ 
nomenological relativistic OP can introduce uncertainties 
in the predictions of the model. The determination of a 
theoretical OP, which fulfills the dispersion relations in 
the whole energy region of interest, would be very useful 
to reduce the theoretical uncertainties. 

In the RGF model the nuclear response is written in 
terms of the single-particle optical-model Green’s func¬ 
tion. This result is obtained retaining only the one-body 
part of the nuclear current. The inclusion of two-body 
MEC would require an extended model based on the two- 
particle Green’s function, whose evaluation represents a 
very hard task. 

The imaginary part of the OP can include pion- 
production processes that should have already been sub¬ 
tracted in the analysis of CCQE data. In this paper we 
have shown that the fully CC-inclusive T2K cross sec¬ 
tions, which include pion production, are clearly under¬ 
estimated by the RGF calculations. Even if we cannot 


disentangle the pion-production contribution that can be 
included in the phenomenological OP, this is not enough 
to reproduce the CC-inclusive T2K data. If we consider 
that the RGF model was developed to describe FSI in the 
inclusive QE scattering and that it is able to give a rea¬ 
sonably good agreement with QE electron and CCQE 
neutrino-scattering data, the present comparison with 
T2K data can be interpreted as an indication that the 
pion-production channel gives only a minor contribution 
to the RGF results. 

The full and explicit inclusion of multinucleon chan¬ 
nels is required to successfully reproduce CC-inclusive 
data. Other models, which explicitly include multinu¬ 
cleon emission channels, obtain a satisfactory agreement 
with the T2K data [55l ]. 
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